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Abstract 


A set of aeroelastic equations describing the 
motion of an arbitrarily mistuned cascade with 
flexible, pretwisted, nonuniform blades is devel- 
oped using an extended Hamilton's principle. The 
derivation of the equations has its basis in the 
geometric nonlinear theory of elasticity in which 
the elongations and shears are negligible compared 
to unity. A general expression for foreshortening 
of a blade is derived and is explicitly used in the 
formulation. The blade aerodynamic loading in the 
subsonic and supersonic flow regimes is obtained 
from two-dimensional, unsteady, cascade theories. 
The aerodynamic, inertial and structural coupling 
between the bending (in two planes) and torsional 
motions of the blade is included. The equations 
are used to investigate the aeroelastic stability 
and to quantify the effect of frequency mistuning 
on flutter in turbofans. Results indicate that a 
moderate amount of intentional mistuning has enough 
potential to alleviate flutter problems in un- 
shrouded, high-aspect-ratio turbofans. 
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Nomenclature 


area centroid and elastic axis 
offset 

unit vectors along x,y,z axes 
unit vectors along x 3 >y 3> z 3 
axes 

shear modulus of elasticity 
nondimensional amplitudes of 
generalized coordinates 
associated with the bending 
modes of the sth 
blade measured out of and in 
the plane of rotation 
nondimensional amplitudes of 
generalized coordinates 
associated with the bending 
modes W- in the rth mode 
of a tuned cascade 


[A],[A r ] 


A J 

a 

a o 

B-^ , B2 
b ’ b R 

c 

E 

[E],[E $>r ] 


cross sectional area of the 
blade 

reference value of A 
aerodynamic matrices 
torsional mode shape 
elastic axis location 
speed of sound 
blade sectional constants 
semi chord and reference 
semichord 
blade chord 

Young's modulus of elasticity 
matrices 

mass and elastic axis offset 
(also base for natural 
logarithm) 


J xx' *zz 


xxo 


[1] 


J 

J 

i 
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bending moment of inertia 
about the major (parallel 
to the x-axis) and minor 
axis through centroid 
reference bending moment of 
inertia 
unity matrix 
%/Tl 

torsional stiffness constant 
reference value for J 
reduced frequency, w Q b/V eff = 

“o b/M eff a o 

polar mass radius of gyration 
about elastic axis 

Ui = £ ♦ kM 

y m m^ J 
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reference polar radius of 
gyration of cross sectional 
mass about elastic axis 
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[M],[M $ ] 
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m o 
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P.P‘ 

[P] 

P J 

[Q] 


mass radii of gyration about 
x and z axes 
polar radii of gyration of 
cross-sectional area about 
elastic axis 
blade length 

lift per unit span, positive 
up (negative z direction) 
aerodynamic coefficients, in 
the rth mode of cascade 
nondimensional lift coeffi- 
cients due to bending and 
torsional motions in the rth 
cascade mode 

nondimensional moment coeffi- 
cients due to bending and 
torsional motion in the rth 
cascade mode 
axial Mach number, V / a rt 

a U 

aerodynamic moment per unit 
span about the elastic axis, 
positive nose up 
relative Mach number, 

/ V a + a " r V a o 

effective relative Mach 
number, V eff /a Q 

number of generalized coordi- 
nates with the bending 
motions out of and in the 
plane of rotation and with 
the torsional motion 
inertial coefficients of the 
sth blade 
inertial matrices 
mass per unit length 
reference mass per unit length 
number of blades in cascade 
an arbitrary point on the 
elastic axis before and 
after deformation 
stiffness matrix, Eq. (16) 
quantity defined by Eq. (4) 
matrix, Eq. (16) 
hub radius 
blade tip radius 
integer (0,1,2, ... N-l) speci- 
fying mode of a tuned cas- 
cade; also blade coordinate 
along elastic axis before 
deformation 
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Wj_ (i = 1,2,...) 
[W],[W] 


xyz 


X 3*3 Z 3 


{ X s ) , {X} 
{Y r ),{Y} 


running blade coordinate along 
elastic axis before and 
after deformation 
position vector of an arbi- 
trary point on the blade 
before and after deforma- 
tion, Eqs. (A6) and (A7) 
position vector of a point on 
elastic axis before and 
after deformation, Eqs. (A6) 
and (A7) 

stiffness coefficients of the 
sth blade 

stiffness matrices 
integer specifying blade, 
s = 0,1,2, N-l; also blade 
spacing 

kinetic energy 

— 2 2 
blade tension, T c = T c /m o n L 

transformation matrix 

time, initial time, final time 

strain energy 

radial foreshortening 

deformations of elastic axis 

in X„, Y a and Z fl 

directions 

axial velocity 

effective relative velocity, 

+ ^ r2 

• cos[90 - s - tan -1 ( V a / « r ) J 

axial extension of elastic 
axis 

work done by aerodynamic 
loading 

beam functions 
model function matrices 
hub-fixed axis system, rotates 
about the Z ft -axis with an 
angular velocity ft 
blade fixed axis system at 
arbitrary point on elastic 
axis 

blade fixed axis system in the 
deformed configuration ob- 
tained by rotating xyz 
column matrices 
column matrices 


j 


a angle of twisting deformation, 

positive when leading edge 
is upward 

a$i (i = 1,2,...) amplitudes of generalized 

coordinates associated with 
torsional modes A- of the 
sth blade 

amplitudes of generalized 
coordinates associated with 
torsional modes l\- in the 
rth mode of tuned cascade 
interblade phase angle in the 
rth mode of tuned cascade, 
2irr/N 

nondimensional eigenvalue, 

(w/(fl 0 )^ 

engineering strain components 
variation of ( ) 
tensor strain components 
damping ratios of sth blade 
blade running coordinate mea- 
sured from hub; n = n/L; 
blade elastic axis position; 

n = (a + 1 ) / 2 

3 2 
reference mass ratio, m /up b D 
o a k 

real part of eigenvalue 
pretwist angle 
fluid and blade material 
density 

nondimensional time, tw 

o 

imaginary part of eigenvalue 
rotational speed 
frequency and reference 
frequency 

angular velocity vector ex- 
pressed in xyz system 
curvature vectors 

derivative a( ) / at or 

3( )/ 3t 

derivative a( )/a r or 
3 ( ) /an 

I. Introduction 

A research program in propulsion system aero- 
elasticity is being conducted at the NASA Lewis 
Research Center. As a part of this general pro- 
gram, an effort was made by the authors in Refs. 1 
and 2, to improve the physical understanding of 
turbofan engine aeroelastic characteristics includ- 
ing blade mistuning (nonidentical blade properties) 
effects. Other published work on mistuning was 
cited in Refs. 1 and 2. 


The mathematical formulation considered in 
Ref. 1 is a 'typical section' model with two 
degrees of freedom, a bending and a pitching about 
the elastic axis. This model was found to be 
sufficient to elicit physical understanding of mis- 
tuning effects and to conduct parametric studies. 
Furthermore, this model was utilized in Ref. 3 to 
show that mistuning has enough potential to sig- 
nificantly raise the flutter speed of an advanced 
fan. This potential may have a very practical 
significance in eliminating the commonly used mid- 
span shrouds in advanced turbofan designs. The 
main purpose of the midspan shrouds is to increase 
the blade natural frequencies and thus to avoid 
aeroelastic instabi 1 ities. However, the shrouds 
have an adverse effect on aerodynamic performance. 

While the typical section model served the 
intended purposes, it is expected to be inadequate 
to obtain accurate flutter boundaries of modern 
technology turbofans or advanced turboprop blades 
which have a considerable amount of pretwist. As a 
consequence of pretwist, bending out of the plane 
of rotation couples elastically with bending in the 
plane of rotation. Consequently, a more realistic 
structural model of a blade is required. 

The purposes of the research summarized in 
this paper are: (1) to develop a more refined and 

realistic structural model of a blade than the 
typical section model and to combine this model 
with the available unsteady cascade aerodynamic 
models; and (2) to study the effects of mistuning 
and other cascade parameters on flutter character- 
istics of an advanced technology fan stage. To the 
best of authors' knowledge, the mentioned effects 
utilizing the purposed model have not been studied 
in the published literature. 

There exist several levels of approximations 
in structural theory used in rotary wing aero- 
elasticity to represent a structural model of a 
blade. The theory of pretwisted blades presented 
in Ref. 4 may be viewed as the first level of 
approximation. The structural theories presented 
in Refs. 5-9 may be viewed as ones with the second 
level of approximation. Since the corresponding 
unsteady cascade aerodynamic theory consistent with 
the second level structural theories is not avail- 
able, the theory of Ref. 4 is used as a first step. 

The governing equations of motion of a 
rotating, pretwisted blade are derived by using an 
extended Hamilton's principle. The derivation is 
characterized by the use of an axial (along the 
beam axis) displacement which includes second- 
degree nonlinear terms defining the axial fore- 
shortening^ of the tension axis due to bending, 
torsion, and the noncoincidence of the elastic axis 
and tension axis. The development of a more gen- 
eral expression for foreshortening than that used 
in Ref. 10 and the explicit use of this general 
expression in the axial displacement field are 
believed to be new. The explicit consideration of 
foreshortening in deriving the equations is appeal- 
ing in that it parallels the corresponding develop- 
ment of equations for a nonrotating blade where 
retention of second-degree terms in the kinetic 
and potential energy expressions leads to linear 
equations (See Ref. 10). A more detailed discus- 
sion on the role of the general expression for 
f oreshortening in deriving the second-degree non- 
linear equations of motion and in accounting for 
the geometric stiffness effects in beam analyses 
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using finite-element methods will be presented in a 
future publication. The disk is assumed to be 
rigid. The unsteady, two-dimensional, cascade, 
aerodynamic loads are calculated by Smith's 
theory 11 in subsonic flow and Adamczyk and 
Goldstein's theory 1 ^ in supersonic flow with a 
subsonic leading edge. The governing equations of 
motion of a mistuned cascade are formulated by 
assuming that the general motion of a blade is a 
linear combination of its motions in all possible 
modes of a tuned cascade. The space variable in 
the coupled integro-partial differential equations 
of motion is eliminated by using a modified 
Galerkin's method. The resulting equations are 
cast as a standard complex eigenvalue problem from 
which the aeroelastic stability is determined. 

A digital computer program is developed to 
form and solve the complex eigenvalue problem. 

This program is written in a modular form such that 
it can be applied to advanced turboprops, heli- 
copter rotors in hover, and wind turbine rotors by 
incorporating the appropriate aerodynamic modules. 
Due to length constraints of the paper, limited 
results are presented only for an advanced turbofan 
blade. However, additional parametric results will 
be presented in a future publication. 

II. Theory 

The components of a bladed disk system have 
complex geometries. The analysis of this system is 
further complicated by blade mistuning. To inves- 
tigate the aeroelastic stability of a cascade with 
mistuning, a mathematical model will be developed 
in this section. 

A. Coordinate Systems 

Several coordinate systems will be employed in 
the derivation of equations of motion;, those which 
are common to both the structural and aerodynamic 
aspects of the derivation for the sth blade are 
shown in Figs. 1 and 2. The axis system 
shown in the figures rotates with a constant angu- 
lar velocity ft about the X n -axis. The Y ft -axis 
coincides with the undeformed elastic axis of the 
blade. The blade principal axes, x and z, of the 
cross section at an arbitrary point on the elastic 
axis are inclined to the X ft and axes by an 
angle s as shown in Fig. 2. The blade elastic 
deformations, u, v, w, and a translate and rotate 
the xyz system to the X 3 y 3 Z 3 system as shown 
in Fig. 1. 

B. Tuned Cascade Model 


If the blades are tuned, the N-bladed cascade 
has N interblade phase angle modes with a con- 
stant phase angle 8 r between adjacent blades. 
This interblade phase angle is restricted by 
Lane's 1 ^ assumption to the N discrete values 
= 2irr/N were r = 0,1,2, ... ,N-1. In each of 
tnese modes all blades have the same amplitude. 

For a tuned cascade, the modes with different 
interblade phase angles are uncoupled. The blade 
deflections are expressed in a traveling wave form 
in terms of a set of generalized coordinates which 
are associated with the nonrotating uncoupled beam 
modes in pure bending and torsion. The number of 
modes retained in the plane of rotation, in the 
plane perpendicular to the plane of rotation, and 
in torsion are M^, Mg, and M tt respectively. 

The blade deflections expressed in traveling wave 
form are 
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An interblade phase angle e r in the range 
0-180° represents a forward traveling wave - a wave 
traveling in the direction of rotation. 

In Eq. (2) the standard nonrotating ortho- 
normal modes for a beam with fixed-free boundary 
conditions are given by 


Wj(n) = cosh(pjn) - cos(pjn) 


(cos p. + cosh p.) 
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A j ( n ) = V? s i n [ ( 2 j - 1) W2 n] 

where the value of pj is obtained from 

cos Pj cosh Pj + 1 = 0 (4) 

In the case of a tuned cascade it is adequate 
to analyze the motion of a cascade in each of the 
interblade phase angle modes separately. Hence, 
the total number of degrees of freedom of a tuned 
cascade is (M^ + Mg + M a ) for each value of B r . 

C. Mistuned Cascade Model 


In the case of an arbitrarily mistuned cas- 
cade, the blades can have different amplitudes, and 
the phase angle between adjacent blades can vary. 
However, the general motion of a blade in a mis- 
tuned cascade can be expressed as a linear combina- 
tion of the motions in all possible interblade 
phase angle modes of the corresponding tuned cas- 
cade. Consequently, the blade deflections can be 
written as 
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For a case with N mistuned blades, Eq. (5) 
can be generalized as 
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The total number of degrees of freedom for 
this general case is N times (M^ + Mg + M a ). 

D. Aerodynamic Model 

The unsteady, two-dimensional, cascade, aero- 
dynamic loads were calculated by using Smith's 
theory in subsonic flow, and Adamczyk and 
Goldstein's^ theory in supersonic flow with a 
subsonic leading edge. In these theories, the air- 
foil thickness, camber and steady state angle of 
attack are neglected, and the flow is assumed to be 
isentropic and i rrotational . At any radial station 
the relative Mach number is a function of the 
inflow conditions and the rotor speed. Most cur- 
rent fan designs have supersonic flow at the tip 
and subsonic flow at the root. As a result, some 
region of the blade span encounters transonic 
flow. Since the above unsteaay theories are not 
valid in the transonic region, the subsonic the- 
ory with M e ff = 0.9 for stations in the range 
0.9 <_ Mpff < 1 and the supersonic theory with 
M e ff = 1.1 for stations in the range 1.0 £ M e ff £ 1.1 
were used. The motion-dependent aerodynamic lift 
and moment were expressed in Ref. 1 in terms of 


nondimensional coefficients. The expressions for 
lift and moment per unit span as a function of the 
same coefficients can. be written as 




The coefficients l hhr , l hor , ..., 1 r are 
calculated for specified values of M e ff, k, s/c, 
£, and a. In the aerodynamic theories used 
herein, the steady state angle of attack is neg- 
lected. In applying these theories to the present 
case, the blade steady state angle of attack is 
also set to zero and the effective relative veloc- 
ity is assumed to be the component of the blade 
relative velocity along the blade chord as indi- 
cated in Fig. 2. 


E. Structural Model 


The structural model of each blade consists of 
a straight, slender, twisted, nonuniform elastic 
beam with a symmetric cross section. The elastic 
axis, the inertia axis, and the tension axis are 
taken to be noncoincident. The effect of warping 
is not explicitly considered. However, a partial 
effect of warping enters into the equations of 
motion as it does in Ref. 4. The blade is assumed 
to be rigid in the direction along the elastic 
axis. Consequently, the axial equation of motion 
is eliminated. The structural model has its basis 
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in the geometric nonlinear theory of elasticity in 
which elongations and shears are negligible com- 
pared to unity and the squares of the derivatives 
of the extensional deformation of the elastic axis 
are negligible compared to the squares of the bend- 
ing slopes. This level of the geometric nonlinear 
theory of elasticity is required to derive a set of 
linear coupled bending-torsion equations of motion. 

F. Equations of Motion 

The equations of motion will be derived by 
using the extended Hamilton's principle in the form 
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The expressions for strain energy, kinetic 
energy, and virtual work of aerodynamic forces can 
be written as 
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The expressions for the strain components, 
position vector, foreshortening, and transformation 
matrix are developed in Appendix A. Substituting 
Eqs. (A9b) , (A9c), and (A14) into Eq. (10), Eq. 

(13) into Eq. (11), Eqs. (8a) and (8b) into Eq. 

(12), the resulting equations into Eq. (9), taking 
the indicated variations, integrating over the 
cross section of the blade wherever necessary, 
integrating by parts over time, neglecting Coriolis 
and rotary inertia terms, and retaining only first- 
degree terms in u, w, and a yields 


The sectional properties in these equations are 
defined in Appendix B. 

By substituting Eqs. (8a) and (8b) into Eqs. 
(14a), (14b), and (14c), nondimensional izing the 
resulting equations, applying Galerkin's method, 
and extending the resultant equations to all the 
blades by using Eq. (6), the equations of motion of 
an arbitrarily mistuned cascade can be simplified as 
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As indicated in Eq. (B2), structural damping 
is added by multiplying the direct stiffness 
coefficients by (1 + 2ichsl)> U + 2 Khs2)> ••• 

(1 + 2ic«s2) ••• 


III. Results and Discussion 


A. Solution 


The aeroelastic stability boundaries are ob- 
tained by solving the standard complex eigenvalue 
problem represented by Eq. (15). The relation 
between the frequency w and y is 

i = iv ^ = p ± i v (17) 

U) 


5. For the mode of interest check whether the 
imaginary part of the eigenvalue is within 
an acceptable tolerance of the reference 
frequency. If not, go back to step (3) 
with a new reference frequency. 

6. For the mode of interest, check whether the 
real part of the eigenvalue is within an 
acceptable tolerance of unity. If not, 
modify the rotational speed and repeat the 
process from step (2). If positive, reduce 
speed; if negative, increase speed. 

7. To find the flutter boundary for another 
value of the axial Mach number go back to 
step (1). 

It should be noted that the above method 
utilizes a rotor speed iteration within the axial 
Mach number iteration loop. It would appear that 
it would be more efficient to place the Mach number 
iteration within the rotor speed iteration loop, 
since the mass and stiffness matrices have to be 
updated each time the rotor speed is changed. 
However, it was found that the time required to 
assemble these matrices is small relative to the 
time required to assemble the aerodynamic matrix 
and to extract the eigenvalues and eigenvectors. 

In addition, the flutter boundary was more 
"distinct" using the method described above. 

The correctness of the program was checked by 
constructing a hypothetical blade model with con- 
stant structural and aerodynamic properties along 
the span. The blade frequencies obtained by using 
the program were checked against known solutions. 
The eigensolutions with aerodynamics included were 
checked against the results obtained by using the 
typical section model program (MISER) developed in 
Refs. 1 and 2. 

For a cascade with nonuniform blades, there 
are no published theoretical results on flutter. 
However the results from ASTROMIC are compared with 
those from MISER by choosing the properties for a 
typical section at different spanwise stations. 
These results will be discussed later. 


Flutter occurs when 7 > 0. 

B. Computer Program and Verification 

A computer program (ASTROMIC) was written to 
assemble and solve the generalized eigenvalue prob- 
lem given in Eq. (15). The program can be used to 
predict the in vacuum natural frequencies of a 
nonuniform rotating beam by setting the aerodynamic 
matrix [A] to zero. In this case the problem is 
reduced to a standard eigensolution. For the 
flutter problem, the aerodynamic matrix is a func- 
tion of the eigenvalues. Hence, an iterative solu- 
tion is required. The iterative technique used 
herein to calculate flutter boundaries (once the 
problem is completely described) is briefly summa- 
rized as follows: 

1. Select an axial Mach number. 

2. Select a rotation speed and assemble the 
stiffness and mass matrices. 

3. Calculate the variation of the reduced 
frequency based on an assumed reference 
frequency a> 0 and relative Mach number 
with span and construct the aerodynamic 
matrix. The initial value for the refer- 
ence frequency is the same as the natural 
frequency of the mode of interest. 

4. Solve the eigenvalue problem. 


C. Aeroelastic Stability of an Advanced Fan 

An advanced unshrouded fan stage (aspect 
ratio = 3.3), representative of a next generation 
fan, was chosen for analysis. A similar stage was 
analyzed in Ref. 3 by using the typical section 
model. The results in Ref. 3 indicated that the 
fan design without shrouds did not meet the flutter 
requirements, but that it may be feasible to use 
mistuning as a passive flutter control. The moti- 
vation for choosing this fan stage for analysis 
herein is to further investigate the use of mis- 
tuning as a passive flutter control. 

The properties of the fan blade are listed in 
Table I and the reference properties are listed in 
Table II. The design point is at a rotor speed of 
4267 rev/min. and an axial Mach number of 0.55. 

At these conditions the tip relative Mach number 
is 1.45. 

The analyses are performed using two modes 
each in the plane of rotation, in the plane per- 
pendicular to the plane of rotation, and in tor- 
sion. Since there are 28 blades in this stage and 
hence 28 interblade phase angle modes, the total 
number of degrees of freedom for an arbitrarily 
mistuned cascade is 168. However, for alternate 
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mistuning, in which every other blade is identical, 
certain symmetry properties are exploited in the 
analysis. Because of symmetry, the 8 r mode 
couples with the (e r + *) mode only. The analy- 
sis for this case consists of 14 separate eigen- 
solutions of 12 degrees of freedom each rather than 
one eigensolution of 168 degrees of freedom for an 
arbitrarily mistuned cascade. 

The blade was analyzed for vibration to gen- 
erate a Campbell diagram which is shown in Fig. 3. 
Only the first and second bending frequencies and 
the first torsion frequency are shown. Since the 
mass ratio for this blade is high, the aerodynamic 
forces result in flutter frequencies which are not 
significantly different from those in a vacuum. 

In Refs. 1-3, a typical section model, which 
was originally used for fixed wing flutter analy- 
sis, was adapted for rotating blades. In the case 
of a fixed wing, the relative velocity of each 
strip along the span is constant; the structural 
properties of the typical section are obtained by 
their respective values at the 3/4 blade span 
station. In the case of a rotating blade, the 
relative velocity of each strip along the blade 
varies. Then, a question arises as to which span- 
wise station should be used as a reference to cal- 
culate the structural and aerodynamic properties of 
the typical section. To answer this question, the 
eigenvalues of the tuned cascade at the design 
condition from both the present beam and typical 
section analyses using various spanwise locations 
as reference are compared in Fig. 4. As expected, 
the cascade is unstable over a range of interblade 
phase angles corresponding to forward traveling 
waves. It should be noted that, whereas there are 
168 eigenvalues, only the predominantly torsional 
eigenvalues are shown. The other modes were found 
to be stable at the design condition. Also, it can 
be seen that the typical section model correspond- 
ing to approximately the 7/10 span reference sta- 
tion gives the best correlation for the unstable 
modes with the nonuniform blade model. This 
finding is in close agreement with the 3/4 span 
reference station usually used for fixed wing aero- 
elastic calculations. This observation is very 
useful in performing preliminary aeroelastic 
analyses. 


By using the iterative procedure defined above 
the flutter boundaries for both the tuned and mis- 
tuned cascades were calculated. The results are 
shown in Fig. 6. For. the tuned cascade it is seen 
that the design point lies well inside in the un- 
stable region. As expected, the axial Mach number 
at flutter monotonical ly decreases with increasing 
rotor speed. For the tuned cascade it should also 
be noted that along the flutter boundary the tip 
relative Mach number is approximately 1.15. 
Consequently, a region near the blade tip may 
experience nonlinearities associated with transonic 
flow which are not accounted for by the unsteady 
aerodynamic theories used herein. The relative tip 
Mach number at the design conditions is 1.45. 

The inclusion of an alternate frequency mis- 
tuning of approximately 7 percent has significantly 
increased the cascade stability by moving the 
flutter boundary to the right. Furthermore, the 
design point is stable for this level of mistuning. 
The mode defining the boundary has maximum partici- 
pation of the even-numbered (or low frequency) 
blades. The tip relative Mach number along the 
boundary of the mistuned cascade is 1.53. Con- 
sequently, the transonic region is considerably 
inboard of the tip and the transonic effects should 
be minimal. 

Figure 6 also illustrates the effect of struc- 
tural damping on flutter of a tuned cascade. A 
structural damping ratio of 0.002 is included in 
each of the torsional modes. As can be seen, the 
damping has a significantly stabilizing effect on 
flutter speed of a tuned cascade. 

IV. Conclusions 

The major conclusions from this investigation 
are summarized as follows; 

1. A general expression for f oreshortening of 
a blade was derived and was explicitly used in 
deriving the first-degree linear equations of 
motion of a rotating blade. 

2. An aeroelastic model and an associated 
computer program for a mistuned cascade with non- 
uniform blades were developed. 


Previous publ ications^ - *^ using the typical 
section model have shown that torsional frequency 
mistuning could have a significant stabilizing 
effect on the cascade. The effect of mistuning is 
further investigated herein with the nonuniform 
blade model. The method used to vary the frequency 
from blade-to-blade was simply to vary the tor- 
sional stiffness, J. For alternate mistuning, the 
torsional stiffness of the odd numbered blades was 
increased by 10 percent over that of a tuned blade 
at each spanwise location. Likewise, the torsional 
stiffnesses of the even numbered blades was de- 
creased by 10 percent from that of a tuned blade. 
The result was a total torsional frequency varia- 
tion of approximately 7 percent. The eigenvalues 
for this mistuned rotor at the design conditions 
are shown in Fig. 5, along with the corresponding 
tuned values. These predominately torsional eigen- 
values have split into high and low frequency 
families corresponding to modes with major partici- 
pation of the odd and even blades, respectively. 

As can be seen, this type of mistuning has stabi- 
lized the cascade to such an extent that it is 
stable at the design point. 


3. For the blade analyzed herein a typical 
section model corresponding to the 7/10 span refer- 
ence station gives the best correlation with the 
nonuniform blade model. 

4. An advanced, unshrouded, high-aspect ratio 
fan was modelled and analyzed with and without mis- 
tuning. The results show that a moderate amount of 
mistuning has enough potential to alleviate flutter 
problems in unshrouded turbofans. 

Appendix A 

Strain Displacement Relations 

A schematic representation of the deformed and 
undeformed elastic axis is shown in Fig. 1. Since 
the pretwist angle K varies with r, with ini- 
tial curvature of the elastic axis is 
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xyz 
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The elastic deformations translate and rotate 
the triad xyz to X 3 y 3 Z 3 . Let the expressions 
for the curvature vector of the deformed elastic 
axis and for the transf ormation matrix between the 
x 3 y 3 * 3 - and xyz systems be 


a) =uje + a) e + to e 

x 3 y 3 Z 3 x 3 x 3 y 3 ^3 z 3 z 3 
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Following the procedure presented in Ref. 8, one 
can develop the second-degree expressions for the 
curvature components and for the transformation 
matrix in terms of u, v, w, and a as 

a) = (u" + aw") sin z + (w" - au" ) cos Z 






+ (u‘ sin £ + w‘ cos £) (u 1 cos £ - w‘ sin £) 
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(A4) 


[T] = 
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(A2) 


- z[(u" + aw") sin z + (w" - au") COS z] 
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Invoking the small strain assumption, the 
engineering strains are related to the components 
of the strain tensor according to: 
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= 2 C 


Y yy e yy’ Y yx € yx’ Y yz € yz 


( A10) 


u' cos Z - w' sin z 
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u* sin z + w' cos c 
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- (u ' sin z + w' cos Z ) 


Z i o 

1 - | ^ ( u ' sin Z + w' cos z) 
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From Fig. 1, the position vector of an arbi- 
trary point on the blade before deformation is 

F 0 = R 0 ( r ) + x * x + z ®z ( A6 ) 


and that of the same point after deformation is 


r i = R 1 ( r i ) + Xe v + 2e , ( A7 ) 

ill x 3 z 3 

Then, the Green's strain tensor based on a 
Lagrangian description is given by 


dr l * dr i “ dr o * dr o = 2 £ dx dr dz Jt £ ij] \ dr \ (A8) 


Substituting Eqs. (A6) and (A7) into Eq. (A8) 
and retaining terms up to second-degree in u, v, 
w, and a, one obtains the following expressions 
for the required components of strain: 


It should be pointed out that in arriving at 
the expressions given in Eq. (A9) several higher- 
order terms have been discarded based either on 
considerations related to small deformations or an 
approximations due to slenderness of the blade. 

It is convenient to eliminate the axial equa- 
tion of motion. This is done by explicitly con- 
sidering the foreshortening due to bending and 
torsion. The expression for foreshortening for 
the present case is derived by making use of the 
equilibrium condition that the integral of the 
longitudinal stress over the cross section must be 
equal to the total tension. Thus, 


// 


yy 


dx dz 


= Ea|v* + i ( u ' 2 + w' 2 ) - e A [(u" + aw") cos z 
- (w " - au") sin £] + k 2 (^- «'!') (All) 
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from which 


v 1 


= v 
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where 
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EA 
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An alternate expression for T c results from 
the kinetic energy and is given in Appendix B. 
Substituting Eqs. (A12) and (A13) into Eq. (A9a) 
and invoking the assumption that the blade is rigid 
(EA » ®) along its elastic axis, one obtains 


'yy 


-z[(u" + aw") sin £ + (w" - au") COS C] 


- ( x - e A )[(u" + aw") cos £ - (w" - au") sin c] 

+ * X 2 - kj) (a' 2 - 2a'C) (A14) 

The other components of strain given by Eqs. 
(A9b) and (A9c) remain unchanged. 

Appendix B 


Sectional Properties and 
Definition of Matrices 
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Note the expression for J is valid only for a 
circular section. For a noncircular section a 
warping correction must be included. 
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TABLE I. - BLADE SECTIONAL PROPERTIES 



n 


0.0199 

0.1017 

0.2372 

0.4083 

0.5917 

0.7628 

0.8983 

0.9801 

EI xx 

0.9719 

0.9665 

0.7827 

0.6072 

0.3067 

0.1633 

0.1149 

0.0839 

EI Z Z 

112.4 

124.8 

152.4 

139.4 

122.7 

119.3 

127.7 

124.1 

GJ 

0.0719 

0.9665 

0.7827 

0.6072 

0.3067 

0.1633 

0.1149 

0.0839 

EA 

0.9974 

1.0225 

0.9711 

0.9358 

0.7641 

0.6482 

0.6038 

0.5541 

ml 

rol 

0.4436 

0.5338 

0.5675 

0.7275 

0.6901 

0.7694 

0.9458 

0.9734 

eb 2 

0 

0 

0 

0 

0 

0 

0 

0 

e 

0 

0 

0 

0 

0 

0 

0 

0 

e A 

0 

0 

0 

0 

0 

0 

0 

0 


0.5859 

0.6095 

0.6265 

0.6722 

0.6972 

0.7461 

0.7995 

0.8226 

X 

0.0940 

0.0925 

0.0854 

0.0767 

0.0603 

0.0478 

0.0415 

0.0370 


1.0103 

1.0516 

1.0816 

1.1616 

1.2060 

1.2913 

1.3840 

1.4242 

L 

m 

0.9974 

1.0225 

0.9711 

0.9358 

0.7641 

0.6482 

0.6038 

0.5541 

b 

1.010 

1.051 

1.082 

1.161 

1.206 

1.291 

1.384 

1.425 


TABLE II. - REFERENCE QUANTITIES 


N = 28 

p m = 4374 kg/m 3 

bR = 0. 0946 m 

p a = 1 kg/m 3 

Rh = 0.3876 m 

a 0 = 340.3 m/sec 

R t = 1.021 m 

A 0 = 0.0034 m 2 

E = 1.23X10 11 N/m 2 

r = 0.5774 

a 

G = 4.744xl0 10 N/m 2 

0 

a> 0 = varies, rad/sec 

I xxo = 9.19x10-8 m 4 

M h = Mg = M a = 2 

J 0 = 3.676xlO“ 7 m 4 

t = tan” 1 [1-552(77 + e H )] 
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Figure 1. - Blade coordinate systems before 
and after deformation. 



Figure 2. - Coordinate systems of blade cross section. 
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Figure 4. - Comparison of eigenvalues from beam model and typical 
section models: w 0 = 1409.6 rad/sec, Q = 446. 8 rad/sec, 

M ax = 0.55. 
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Figure 5. - Comparison of eigenvalues of tuned and mistuned 
cascade: oj q s 1391.1 rad/sec, Q = 446.8 rad/sec, M ax = 0.55. 



Figure 6. - Comparison of tuned and mistuned flutter boundaries. 
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